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Abstract. In this paper, inverse spectral problems for Sturm-Liouville operators on a tree 
(a graph without cycles) are studied. We show that if the potential on an edge is known a priori, 
then 6 — 1 spectral sets uniquely determine the potential functions on a tree with 6 external 
edges. Constructive solutions, based on the method of spectral mappings, are provided for the 
considered inverse problems. 
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1. Introduction. 

This paper concerns the theory of inverse spectral problems for Sturm-Liouville operators 
on geometrical graphs. Inverse problems consist in recovering differential operators from their 
spectral characteristics. Differential operators on graphs (quantum graphs) have applications in 
various helds of science and engineering (mechanics, chemistry, electronics, nanoscale technology 
and others) and attract a considerable attention of mathematicians in recent years. There is 
an extensive literature devoted to differential operators on graphs and their applications, we 
mention only some research papers and surveys [D El El m El El [ig. 

There are different kinds of inverse problems studied for quantum graphs, one of them is to 
recover the coefficients of the operator while some information is known a priori. This paper is 
focused on the reconstruction of the potential of the Sturm-Liouville operator on a tree (a graph 
without cycles) with a prescribed structure and standard matching conditions in the vertices. 
V.A. Yurko [HE] studied such inverse problems on trees by the Weyl vector, the system of 
spectra and the spectral data. These problems are generalizations of the well-studied inverse 
problems for Sturm-Liouville operators on a finite interval (see monographs piininiE] and 
references therein). By the method of spectral mappings [l2l|T3], V.A. Yurko proved uniqueness 
theorems and developed a constructive algorithm for solution of inverse problems on trees. 

In this paper, we formulate and solve partial inverse problems for the Sturm-Liouville op¬ 
erator on the tree. We suppose that the Sturm-Liouville potential is known on the part of 
the graph and show that we need less data to recover the potential on the remaining part. 
We know the only work m in this direction, where the potential is known on a half of one 
edge and completely on the other edges of the star-shaped graph, and the author solves the 
Hochstadt-Lieberman-type problem [ig by a part of the spectrum. 

In this paper, we assume that the potential is known on one edge of a tree, then reconstruct 
the potential on the remaining part by the system of spectra or the Weyl functions. By 
developing the ideas of V.A. Yurko [ZIE], we show that one needs one less spectral set or 
one less Weyl function for the solution of the partial inverse problem. We consider separately 
the cases of boundary and internal edges, present constructive solutions and corresponding 
uniqueness theorems for both of them. 

The results of this paper can be generalized to the case, when the potential is known on 
several edges. However, in this case the number of given spectra, sufficient to recover the 
potential on the whole graph, depends not only on the number of these edges, but also on their 
location (see the example in Section 5). We note that the method of spectral mappings works 
also for graphs with cycles (see [IS]), so one can generalize our results in this direction. 
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The paper is organized as follows. In Section 2, we introdnce the notation and briefly 
describe the solution of inverse problems on trees by V.A. Yurko EE]. In Section 3, we 
formulate our main results and outline their constructive solutions. Section 4 contains proofs 
of the technical lemmas from Section 3. In Section 5 we illustrate our method by an example. 

2. Inverse problems on a tree 

In this section, we introduce the notation and provide the main results of V.A. Yurko on 
the inverse problems on trees (see works EE] for more details). 

Consider a compact tree G with the vertices V = and edges E = For 

each vertex n G V, we denote the set of edges associated with v by and call the size of E^ 
the degree of v. Assume that the tree G does not contain vertices of degree 2. The vertices of 
degree 1 are called boundary vertices. Denote the set of boundary vertices of the graph G by 
dG. For the sake of convenience, let each boundary vertex Vi be an end of the edge e*, such 
edges are called boundary edges. All other vertices and edges are called internal. Let the vertex 
Vr G dG be the root of the tree. 

Each edge Cj G E is viewed as a segment [0,Tj] and is parametrized by the parameter 
Xj G [0, Tj]. The value Xj = 0 correspond to one of the end vertices of the edge ej, and xj = Tj 
corresponds to another one. For a boundary edge, the end Xj = 0 corresponds to the boundary 
vertex vj. 

A function on the tree G can be represented as a vector function y = [yjj^^i, where yj = 
yj{xj), Xj G [0,Tj], j = l,m. Let Cj = [vi,Vk\, i.e. the vertex Vi corresponds to the end Xj = 0 
and the vertex Vk corresponds to Xj = Tj. Introduce the following notation 

Vji'Vi) = 2/i(0)> Vji'Vk) = VjiTj), 
y'j (d) = y'j (0), y'j {vk) = -y'j {Tj). 

If Vi G dG, we omit the index of the edge and write y{vi) and y'{vi). 

Consider the Sturm-Liouville equation on G\ 

- y'j + qj{xj)yj = Xyj, Xj G [0, Tj], j = Cm. (1) 

where A is the spectral parameter, qj G L[0,Tj]. We call the function q = [qj]'jLi the potential 
on the graph G. The functions yj, y'j are absolutely continuous on the segments [0,Tj] and 
satisfy the standard matching conditions in the internal vertices v G V\dG: 

\yj{v) = yk{v), ej,ek G (continuity condition), 
y'ji'^) = 0; (Kirchhoff’s condition). 

Let Lq and L^, & dG, be the boundary value problem for system ([T|) with the matching 

conditions (E) and the following conditions in the boundary vertices: 

Lq-. y{vi) = t), Vi^dG, (3) 

Lk:y'{vk) = 0, y{vi) = 0, ViedG\{vk}. ( 4 ) 

It is well-known, that the problems Lk have discrete spectra, which are the countable sets of 
eigenvalues A^. = {Afcs}“ i, A: = 0 or G dG. 

Fix a boundary vertex Vk G dG. Let = ['ipkj]'jLi, = '^’kjixj. A), be the solution of the 
system o. satisfying the matching conditions (E) and the boundary conditions 

V^fcfc(0, A) = 1, '0fci(O, A) = 0, Vj G aG\{nfc}. 

Denote Mfc(A) = A). The functions and are called the Weyl solution and the Weyl 

function of (E) with respect to the boundary vertex Vk, respectively. The notion of the Weyl 
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function for the tree generalizes the notion of the Weyl function (m-function) for the classical 
Sturm-Liouville operator on a hnite interval mm If the tree G consists of only one edge, 
then Mfc(A) coincide with the classical Weyl function. 

Consider the following inverse problems. 

Inverse Problem 1. Given the spectra Aq, Vk G 9G'\{nr}, construct the potential q on the 
tree G. 

Inverse Problem 2. Given the Weyl functions Mk{X), Vk G 9G'\{nr}, construct the potential 
q on the tree G. 

Note that if the number of boundary vertices is b, then one needs b spectra or 6 — 1 Weyl 
functions to recover the potential. We do not require the data associated with the root Vr- 

There is a close relation between Inverse Problems [T] and [21 The Weyl functions can be 
represented in the form 

A'4(A) = -h^, vtedG, (6) 

where Afc(A) are characteristic functions of the boundary value problems L^. If the eigenvalues 
Afc are known, one can constrnct characteristic fnnctions as infinite products by Hadamard’s 
theorem. Thus, with the system of spectra, one can obtain the Weyl fnnctions and rednce 
Inverse Problem [1] to Inverse Problem [21 

V.A. Ynrko has proved, that Inverse Problems [T] and [21 are nniqnely solvable, and provided 
a constrnctive algorithm for the solntion by the method of spectral mappings [12]. In the 
remaining of this section, we shall briefly describe his algorithm. Let the Weyl fnnctions Mfc(A), 
Vk G 9G\{nr} be given. Consider the following anxiliary problem. 

Problem IP(k). Given Mk{\), constrnct the potential qk{xk) on the edge Ck- 

Note that this problem is not eqnivalent to the inverse problem on the hnite interval, since 
the Weyl fnnction Mfc(A) contains information from the whole graph. However, it can be solved 
nniqnely by the method of spectral mappings, and the potential on the bonndary edges can be 
recovered. Then V.A. Ynrko used so-called /r-procedure to recover the potential on the internal 
edges. We reformnlate these ideas in the form, which is more convenient for us in the fnture. 

Theorem 1. Let v he an internal vertex, connected with the set of boundary vertices V C 
9G'\{nr} only one other vertex. Suppose the potentials qk on the edges Ck are known for 

all Vk G V', as well as a Weyl function Mk{X) for at least one vertex from the set V. Denote 
G' the graph by removing the vertices Vk G V together with the corresponding edges Ck from 
the graph G. Then the Weyl function for the graph G' with respect the the vertex v can be 
determined from the given information. 

Applying Theorem[T], one can cut the boundary edges off, nntil the potential will be recovered 
on the whole graph. 

3. Partial inverse problems 

In this section, the main resnlts of the paper are formulated. We assume that the potential 
is known on one edge of the tree and formulate partial inverse problems. We consider separately 
the cases of boundary and internal edge. The first one appears to be trivial, for the second one 
we describe the procedure of the constructive solution. For the convenience of the reader, the 
proofs of the technical lemmas are provided in Section 4. 

Inverse Problem 3. Let ej be a boundary edge (f ^ r). Given the potential qj on the edge 
Cf and the spectra Aq, A^, Vk G dG\{vf,Vr]. Construct the potential q on the tree G. 
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The solution of Inverse Problem 3 is a slight modihcation of the method described in 
Section 2 . From Aq, A^, Vk G dG\{vf,Vr}, we easily construct the potentials qk for 
Vk G dG\{vf,Vr}. The potential qj is known, so we can apply Theorem [T] iteratively and 
recover the potential on G. 

Now let 6 / be an internal edge. If this edge is removed, the graph splits into two parts, call 
them Pi and P 2 - Let dPi and dP 2 be the sets of boundary vertices of Pi and P 2 , respectively. 
Fix two arbitrary vertices Vri G dPi and Vr2 G dP2. 

Inverse Problem 4. Given the potential qf on the internal edge Cf, the spectra Aq, A^ Vk G 
9G\{nri, nr 2 }- Construct the potential q on the tree G. 

Solution of Inverse Problem [H For simplicity, we assume that the ends of the edge Cf 
have degree 3. The general case requires minor modihcations. If one splits each of the ends of 
6 / into three vertices, the tree splits into hve subtrees Gi, i = 1,5, such that Vri G G 2 , Vr 2 G G 5 , 
and G 3 contains the only edge e/ (see £g. [I]). Let vi and are arbitrary boundary vertices of 
the trees Gi and G 4 (different from the ends of e/), Vri = V 2 , Vr 2 = v^, Cf = [^ 3 ,^ 6 ]- 



Figure 1 


Step 1. Construct the characteristic functions Afc(A) by the given spectra A^., k = D and 
Vk G dG\{v2,V5}. Find Mk{X) by formula ([5]). 

Step 2. Consider trees Gi and G 4 . Recover the potential q on the edges of Gi and G 4 , 
using the solutions of the problems IP(k) for Vk G (?Gi\{n 3 } and Vk G 9G4\{n6}, and them 
applying Theorem [T] iteratively. 

Step 3. Introduce the characteristic functions of the boundary value problems for the 
Sturm-Liouville equations ([T]) on the graphs G 1 -G 5 with the standard matching conditions ([ 2 ]) 
in internal vertices and the following conditions in the boundary vertices: 


graph Gi 


Af'^(A): y{vk) = 0, Vk e dGi, 

Af^(A): |/'(ni) = 0, y{vk) = 0, VkedGi\{vi}, 

Af^(A): y’{v3) = 0, y{vk) = 0, Vk e ^^^{^ 3 }, 

Af^(A): y'{vi) = 0, y'{vii) = 0, y{vk) = 0, Vk E dGi\{vi,V 3 }. 


graph G 2 
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graph Gs 


graph G 4 


graph G 5 


Af°(A) 

Vivs) 

= 0 , 

yive) = 0, 


A3™(A) 

y'ivs) 

= 0 , 

o' 


Af"(A) 

yivs) 

= 0 , 

y'ive) = 0, 


A3™(A) 

y'ivs) 

= 0 , 

y'{ve) = 0. 


Af°(A) 

y{vk) 

= 0 , 

Vk G dG^, 


A3™(A) 

y'{vi) 

= 0 , 

CD 

Vk G 9G4\{n4}, 

A5>«(A) 

y'iv^) 

= 0 , 

0 " 

Vk G c?G4\{n6}5 

A3™(A) 

y'ivi) 

= 0 , 

CO 

^0 

y{vk) = 0, Vk 

Af(A): 

yi^k) = 

0 , 

Vk G OGq, 


Ar(A): 

y'ive) = 

0 , 

CD 

Vk G 5G5\{n6}. 


Lemma 1. The following relation holds 






Ar(A) -Af(A) 

0 

0 

0 

0 


0 Af(A) 

-1 

0 

0 

0 

Ao(A) = 

Ar(A) A^(A) 

0 

-1 

0 

0 

0 0 

A? 

CO 

< 

Af^(A) 

-Ar(A) 

0 


0 0 

0 

0 

Af^(A) 

-Af(A) 


0 0 

A3^^(A) 

A3^^(A) 

Af^(A) 

A^(A) 


( 6 ) 


If one changes Af'^(A) to A]^^(A) and Af'^(A) to A]^^(A), he obtains the determinant equal 
to Ai(A). Similarly, if one changes Af^(A) to A^^(A) and A|^^(A) to A^^(A), he gets A 4 (A). 

Step 4. Note that the functions Ao(A), Ai(A), A 4 (A) are known from Step 1. Since we 
know the potential on the graphs Gi, G 4 (from Step 2 ) and G 3 (given a priori), we can easily 
construct the characteristic functions for these graphs. Consider the relation ([S]) and similar 
relations for Ai(A) and A 4 (A) as a system of equations with respect to Af'(A), A^(A), Af’(A) 
and ('^) in the following form 


r Oil A^ A^ + ai2 A^ A^ + ai^A^ A^ + ai4A^ A^ = Aq, 

< 021A 2 Af* + a22A^A^ + a2sA2 A^ + a24A^ A^ = Ai, (7) 

[ (^31A 2 A^ + 032A^A^ + 033A^ A^ + as^A^A^ = A4, 


where Oij = aij{X), i = 1,3, j = 1,4, are known coefficients. 

Step 5. Multiply the hrst equation of o by Ai and subtract the second equations, 
multiplyed by Aq. Apply the similar trasform to the hrst and the third equations. Then we 
obtain the system 


611 Af Af + 6 i 2 A^Af + 6 i 3 Af Af + feuA^Af = 0, 

621 Af* A^ + 622 A^ Af + 623Af’ A^ + 624A^ A^ = 0 , 


where 

bii = oiiAi — a2iAo, b2i = aijA4 — osiAo, i = 1,4. 
Divide both equations by A^Af*. 

bii + &i2^2 + bi^M^ + 6j4M2M5 = 0, i = 1, 2, 


where 


M 2 (A) 


Ar(A) 

Af(A)’ 


Ms (A) 


Ar(A) 

Af(A) 


( 8 ) 

(9) 
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are (up to the sign) the Weyl functions for the subtrees G 2 and G 5 associated with the vertices 
V 3 and vq, respectively. 

Step 6. From the system (E]) we easily derive 


Ms 


bil + 

bi3 + &i4^2 


i = 1,2. 


Hence 

(611 + &12d^2)(&23 + &24d^2) = (^21 + ^22d^2)(&13 + &14d^2)- 
Finally, we obtain the quadratic equation with respect to M 2 (A): 

A{X)M^{X) + H(A)M 2 (A) + C'(A) = 0, 

with analytic coefficients ^(A), B{X), C(A): 

A = 612 &24 — ^22^14, 

B = 611624 + 612623 — 621614 — 622613 , 

C = 611623 — 621613 . 


( 10 ) 


( 11 ) 


Step 7. Consider the Sturm-Liouville equation ([T]) on the tree G with the potential g = 0 . 
Implement Steps 1-6 for this case and obtain the quadratic equation 

Ao{X)Ml^{X) + Ho(A)M2o(A) + C'o(A) = 0 , ( 12 ) 

analogous to (ITUD . Denote p = \/A, Rep > 0 , := {p: Rep > 0 , |Imp| < 6}, 6 > 0 , 

[ 1 ] = l+ 0 (p“^). Let /(p^) be an analytic function and £ > 0 . Denote Ze(/) := {p: |/(p^)| > £}. 

Lemma 2. The following asymptotic relations hold 

H(A) = Alo(A)[l], B{X) = Ho(A)[l], C'(A) = C'o(A)[l], p G n Z,{AoBoGo), |p| ^ cx). 

Consequently, D(A) = Do(A)[l] for p E Ss H Z^^Dq), |p| —)■ 00, where D{X) and Do(A) are 
discriminants of equations flTOl) and ([12]), respectively. 

Lemma 3. Ao(A) ^ 0 , Do(A) ^ 0 . 

It follows from Lemmas [ 2 ] and [ 3 l that the quadratic equation flTOl) does not degenerate for 
p E Ss A Z^^AqDq), and two roots of ffTOj) are different by asymptotics as |p| —)■ 00. One 
can easily find an asymptotic representation of M2 (A) for any particular graph and choose the 
correct root of fITU]) on some region of Ss for sufficiently large |p|. Then the function M2 (A) can 
be constructed for all A G C except its singularities by analytic continuation. Similarly one can 
find M5(A). 

Step 8. Consider the tree G2 with the root V2- Solve problems IP(k) by Mfc(A), Vk G 
dG2\{v2, ue}, and by M2(A) for U3, obtain the potential on the boundary edges except 62. Then 
apply the cutting of boundary edges by Theorem [ 1 ] and recover the potential q on G2. The 
subtree G5 can be treated similarly. 

Thus, we recovered the potential q on the whole graph G. In parallel, we have proved the 
following uniqueness theorem. 

Theorem 2. Let the potential qf on the edge eg (f ^ r) be known. 

(i) If ef in a boundary edge, the spectra Aq, A^, vu E dG\{vf,Vr}, uniquely determine the 
potential q on the whole graph G. 

(a) If Cf is an internal edge, the spectra Aq, A^ Vk E 'i^r2} uniquely determine the 

potential q on the whole graph G. 
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Using the described methods with some technical modifications, one can solve partial inverse 
problems by Weyl functions. 

Inverse Problem 5. Let e/ be a boundary edge (f ^ r). Given the potential qj on the edge 
6/ and the Weyl functions Vk G dG\{vf,Vr}. Gonstruct the potential q on the tree G. 

Inverse Problem 6. Given the potential qj on the internal edge Cf, the Weyl functions Mk{X) 
Vk G 9G\{nri, nr 2 }- Gonstruct the potential q on the tree G. 

Thus, if the number of boundary edges is b and the potential is known on one edge (boundary 
or internal), b — 2 Weyl functions are required to construct q on the whole graph. 

4. Proofs 

4.1. Proof of Lemma [T]. Consider the Sturm-Liouville equation ([T]) on the tree G. Let 
Cj{xj, A) and Sj{xj, A) be solutions of ([I]) on the edge Cj under initial conditions 

C,(0, A) = ^'(0, A) = 1, C'(0, A) = ^,(0, A) = 0. 

Any solution y = [yj]fLi of the equation ([T]) on G admits the following representation 

yj{xj,X) = M^{X)Cj{xj,X) + M]{X)Sj{xj,X), j = l,m, xj G [0,Tj]. (13) 

Let BC be some fixed boundary conditions in the vertices v G dG of the form y{v) = 0 or 
y'{v) = 0 (for instance, we consider conditions ([3]) for the problem L and dl]) for the problem 
Lk). Denote by L the boundary value problem for the Sturm-Liouville equation ([1]) with the 
standard matching conditions ([2]) and the boundary conditions BC. If ?/ is a solution of a 
boundary value problem L, substitute d into ([2]) and BC, and obtain a linear algebraic 
system with respect to M°(A), Mj{X). It is easy to check that the determinant of this system 
is a characteristic function A(A) of the boundary value problem L, i.e. zeros of A(A) coincide 
with the eigenvalues of L. 

Example 1. Consider the problem Lq for the star-type graph for m = 3. Then boundary 
conditions ([3]) yield M°(A) = M^^X) = Mg (A) = 0. Consequently, from (|2]) we obtain the 
system with respect to M°(A), j = 1, 2, 3, with the determinant 


Ao(A) 


S,{T,,X) -^2(T2,A) 0 

0 S2(T2,X) -53(T3,A) 

S[{T,,X) A'(T2,A) ^'(T3,A) 


In the general case, the following assertion is valid. 


Lemma 4. Let w ^ V and the degree of w be equal n. Splitting the vertex w, we split C into 
n subtrees Ci, i = l,n. For each i = l,n, let Af (A) and Af^(A) be characteristic functions for 
boundary value problems for equation o on tree Ci with matching conditions (0. boundary 
conditions BC for v G dC fl dCi and the Dirichlet condition y{u) = 0 for Af (A) and the 
Neumann condition y'{u) = 0 for AA^(A). Then the characteristic function A(A) for C with the 
conditions ([2]) and BC admits the following representation: 


A(A) 


Af(A) 

-A^(A) 

0 

0 

0 

Af(A) 

-Af(A) 

0 

0 

0 

0 

■ - A?(A) 

Af(A) 

A^(A) 

A^(A) 

A^(A) 


(14) 
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Indeed, if we write the determinant for A (A) and analyze the participation of the edges 
of Gi in this determinant, we can easily see that A(A) = Af(A)Di(A) + Af^(A)Ei(A), where 
the fnnctions -Dj(A) and Bi(A) do not depend on the subtree G*. Thus we can consider the 
simplest case of the star-type graph, when each Gi contains only one edge, and then change 
the multipliers, corresponding to subgraphs Gi, to Af’(A) and Af^(A). Thus we directly obtain 
fimi from the formula for the star-type graph. 

Lemma [1] follows from Lemma M for the graph in the fig. [TJ Alternatively, one can derive 
(p|) from changing characteristic functions for one-edge subtrees by general characteristic 
function. 

4.2. Proof of Lemma [21 Together with L consider the boundary value problem for 
equation ([T]) with q = 0, the matching conditions ([ 2 ]) and the boundary conditions BG. If some 
symbol 7 denotes the object related to L, we denote by the symbol 7 ° the similar object related 
to In particular, A°(A) is the characteristic function of Let the symbol P{p) stand for 
different polynomials of sin pTj and cos pTj, j = l,m. 

Lemma 5. The characteristic function A(A) has the following asymptotic behavior: 

A(X) = A»(A) + 0(p-‘‘) = ^ + 0(p-‘‘), p e Si, IpI ^ 

pa-L 

where P{p) ^ 0 and d = m — i — n, where m is the number of the edges, i is the number of 
internal vertices and n is the number of boundary vertices with the Neumann boundary condition 
y'{v) = 0. 

Proof. The claim of the lemma immediately follows from the standard asymptotic formulas 
Gj{xj, A) = cospxj + 0{p~^), Gj{xj, A) = —psinpxj + 0(1), 

Sj{xj, A) = + 0(p“^), S'Jx, A) = cospXj, p e Ss, |p| -)■ 00 , 

p 

and the construction of A (A). The relation P{p) ^ 0 follows from the regularity of the standard 
matching conditions. □ 

Applying Lemma |5] to the characteristic functions, dehned on Step 3 of the algorithm, we 
derive asymptotic representations for the coefficients c = Oij, bij, A, B, G in the following form: 

c(A) = c°(A) -I- 0{p '^) = + 0{p p G Ss, IpI —)■ 00 , 

pa -1 

where d stands for different integers. This relation yields Lemma [2l 

4.3. Proof of Lemma [Si In this subsection, we consider only the problem with q = 0, 
so we omit the index 0 for brevity. For simplicity, let Tf = 1. Taking into account, that 

A 3 ™ = -psmp 

and doing some algebra with the expressions (E]), dH}, ffTT]) . we derive 

A(A) = -Fi(A)F4(A)A„(A)2!AAf°(A)Af(A)x(A), (15) 

P 

B(A) = -F,(A)F4(A)5^A„(A) |Af (A)n(A) + Af (A)^«A) - Af“(A)A 3 "(A)^x(A)} . 

(16) 
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( 17 ) 


C(A) = f’,(A)fi(A)A„(A)^Af (A) |n(A) + 5!^5(A)|, 
where 

F,(A) = Af^(A)Ar (A) - Af^(A)Ar (A), ^ = 1,4, 

n(A) = 2Af^(A)A?^(A)Af^(A), 

x(A) and ^(A) are characteristic functions of the graphs Gi U G 2 U G 3 and Gi U G 2 U G 3 U G 4 , 
respectively. Here we mean that the copies of the vertex ^3 (and Ug in the second graph) are 
joined into one vertex with the standard matching conditions () 2 ]). 

Lemma 6. Let Vi and V 2 he two fixed vertices from dG. Denote by A^^{X) 

and A'^^(A) the characteristic functions for equation ([T]) on the tree G with the matching 
conditions o, with the following boundary conditions: 

A^^(A): y{vi) = y{v2) = 0, 

A^^(A): y{vi) =y'{v 2 ) = 0 , 

A^^(A): y'{v,) = yiv 2 ) = 0 , 

A^^(A): y'{v,) = y'iv 2 )= 0 , 

and with the conditions BG in the vertices v G clG\{ni,n 2 }. Then 

ADD(a)a^^(A) - A^^(A)A^^(A) ^ 0 . (18) 

Proof. We shall divide the proof into the following steps. 1. Let the tree G consists of the only 
edge [^ 1 ,^ 2 ]. Then one can check the relation ffTSjl by direct calculation. 

2. Let the vertices Vi and V 2 be connected by edges with the same vertex n, and let there 
also be subtrees Gj, z = l,n, from the vertex v (see £g. H]). Denote by Af’(A) and Af^(A) the 
characteristic functions for Gi with the matching conditions ([2]), the boundary conditions BG 
and y{v) = 0 for Af(X) and y'(v) = 0 for Af^(A). According to LemmalU the following relation 
holds 

A'^'^(A) = —— ^ ^ —^-^A^(A) H—(sinpTi cospT 2 + cospTi sinpT 2 )A^(A), 

P 

where 

A"(A) = nAf(A). ApA) = A"(A)X:^. 

i=l i=l 

Using similar representations for A^'^(A), A^^(A) and A'^^(A), we derive 
ADD(a)a^^(A) - A^^(A)A^^(A) = - {A^{X)Y ^ 0 . 

3. Now let the vertices vi and V 2 be connected by the edges with ^3 and 1 ) 4 , respectively. 
Let the tree G splits by the vertices ^3 and into the subtrees Gj, i = l,ni, connected 
with Us, the subtrees Gj, j = l,'n 2 , connected with ^ 4 , the subtree Gq, including the both 
vertices ^3 and ^ 4 , and the edges Ci, 62 (see £g. [3]). Denote by Af’(A), A^{X), i = l,ni, 
and by Aj^(A), A^(A), j = 1 ,^, 2 , the characteristic functions for the subtrees G* with the 
Dirichlet or Neumann boundary condition in ^3 and for the subtrees Gj with the Dirichlet or 
Neumann boundary condition in ^ 4 , respectively. Let A^'^(A), A^'^(A), A^'^(A) and Aq^{X) 
be characteristic functions for the subtree Gq with the following boundary conditions 

yivs) = y{v4) = 0, 
y{v 3 ) = y'{vi) = 0 , 

A^^(A): y'{v3) = y{v^) = 0, 

Ao^^(A): yfv3) = y'iv,) = 0 , 
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Figure 2 



Figure 3 


and the conditions BC m. other boundary vertices. Denote the functions 


i=l 


A?(A) = nAfW. 

i=i 


Af(A) = AyA)5^ 

i=l 


Af(A) 

Af(A)' 


Af (A) = A?(A) ^ 


Af(A) 

Af(A)' 


A*-''-(A) = Af°(A)Af (A)Af (A) + A™(A)An(A)Af (A) 1 

+Ar(A)Af (A)An(A) + AS'«(A)Aj'(A)An(A), 
A"'f(A)= A°°(A)A?(A)Af(A) + AD''(A)An(A)An(A), I 

A™(A) = Ar (A)Af (A)An(A) + A™(A)An(A)An(A), 
Ann(A) = A^ (A)An(A)An(A). J 

In view of Lemma IU the following relation holds 


( 19 ) 


sinpTisinpTa cos pTi sin pTs ni^ 


A""(A) = 




-A^^(A) + 


P 


A''^(A) 
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sin pTi cos pT2 ^ f w , rp T^Ann^\\ 

H-A (A) + cospTi cospT 2 A (A). 

P 

Together with the similar relations for A'^'^(A), A'^'^(A) and A^^(A), it yields 

A^^(A)A^^(A) - A^^(A)A^^(A) = A™(a)A^^(A) - An^(A)A^n(A) 

Taking flT^ into account, we obtain 

A™(a)A^^(A) - = (A^^(A)A^'^(A) - A^'^(A)A^^(A)) (An(A)An(A))'. 

By virtue of Lemma 0, Ap(A) ^ 0 , z = 1 , 2 . Therefore the relation flTSD holds for the tree G if 
and only if it holds for the subtree Gq. By induction, the claim of the lemma is valid for any 
tree G. □ 

By virtue of Lemmas 0 E] and flT^ . A(A) ^ 0. It follows from flT^ . flTB]) . flTTD . that 
D{X) = B^{X) - AA{X)G{X) 

= F?(A)F 4 ^(A)^A„(A) |Af(A)n(A) + Af(A)^«A) + Af'’(A)Aj'(A)5!!^x(A)|'. 

P I P P ) 

Note that the expression in the bracket above equals to 

Af (A)n(A) + 2 !^a„(a). 

Similarly to Lemma El the following asymptotic formulas can be obtained: 

Af (A)n(A) = C'ir“^exp(r(T - 1))[1], ^^Ao(A) = C' 2 r“'? exp(r(T + 1))[1], 

m 

where p = ir, r ^ +oo, T = ^ Tj, Gi, G^-, p and q are some constants. Clearly, the second 

i=i 

term grows faster than the hrst one. Therefore Ao(A) ^ 0 implies D{X) ^ 0. The proof of 
Lemma [3] is hnished. 

Using Lemma El one can also check, that B(X) and D(X) have the same power of p in the 
denominator, so the roots of f[T 0 |) have different asymptotic behavior. 

5. Example 

In this section, we provide the solution of Inverse Problem E] for the example of the graph 
in the £g. El For simplicity, let Tj = 1, j = 1, 5. Let X 3 = 0 corresponds to the vertex ^3 and 
X 3 = 1 corresponds to Vq. For the boundary edges, Xj = 0 correspond to the boundary vertices. 
The matching conditions ([2D take the form 

V3: z/i(l) = z/2(l) = 2/3(0), 2/'i(l)+ 2/^(1)-2/^(0) = 0 , 4 . 

^^ 6 : 2/3(1) = 2/4(1) = 2/5(1), 2/3(1)+2/4(1)+ 2/^1) = 0 . 

For this example, each subtree Gi consists of only one edge e*, i = 1,5. Let us know 
the spectra Aq, Ai, A 4 and the potential q^. Using the given spectra, one can easily hnd the 
characteristic functions Ao(A), Ai(A), A4(A) and the Weyl functions Mi (A), M^^X). Solving 
problems IP(1) and IP (4), recover qi and q^. 
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Figure 4. Example 


Va 


V5 


Consider the boundary value problem L. Represent the solution y in the form flT^ and 
substitute it into ([2]) and ([3]). From ([3]), one gets M°(A) = (A) = M2{X) = M^{X) = 0. 

Then matching conditions fl20|) yield the system 


/5i 

-52 

0 

0 

0 

0 \ 


(Mix 

0 

52 

-1 

0 

0 

0 


Ml 

51 

51 

0 

-1 

0 

0 


Ml 

0 

0 

C'3 

53 

-54 

0 


Ml 

0 

0 

0 

0 

54 

-55 


Ml 

VO 

0 


51 

51 

S'J 


WJ 


( 21 ) 


Here we omit arguments (1,A) and (A) for brevity. The characteristic function Ao(A) equals 
the determinant of m- Since we know qi, and ^ 4 , we can solve ([I]) and obtain the func¬ 
tions Sj{xj,X) and Cj(xj,X) for j = 1,3,4. Therefore the determinant admits the following 
representation 

Aq = aniS2S5 -f auS^S^ + ai3S'25'5 -1- a^S^S'^, 


where 


an = R; 


^3 -^4 

S' R' 


+ 5i 


ai3 = (51^3 + ^3)54, 


^3 -54 

C' 51 

ai4 = 5 i5354. 


ai2 — 5i 


53 -54 


51 


51 


If one change 5i to Ci or 54 to C4, he obtains analogous relations for Ai(A) and A4(A), 
respectively. Thus we arrive at the system ([7]). 

Let g = 0 on G. Then 


Cj(xj,X) = cos pxj, Sj(xj,X) = 


sm pXj 
P 


Q sin3p Q Q sin2psinp „ sin p 

an = —, On = a^g =--, 044 = — 


P 

,0_0 


®21 ~ ®31 ~ ^Pi ®22 ~ ®33 ~ 


Q sin 2 pcosp 




P 


g Q cos 2 psinp Q Q cos p sin p 

— , a24 = a34 = — . 


®23 ~ ®32 ~ 


p 


P" 
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—9sin5p + 13sin3p + 6sinp —9 cosSp + 7cos3p + 2 cosp 

— IT ~ 5 — ^ 4 . — 


16p^ 


16p2 


Using dH]), we obtain 


—3sin6p — 2sin4p + 13sin2p —3 cos bp + 6 cos 4/? + 3 cos 2/? — 6 

"11 = "21 = ~ ~ 


16p3 


_ LU _ 
, C]_2 — ^23 — 




^g 3 cos 6 p — 10 cos4p + 13 cos2p — 6 ^g ^g —3sin6p + 12sin4p — 15sin2p 

Dio = Don = 


32p5 


, „ _: fv -}V - 

^13 ~ *^22 — 22 p 4 ’ — "24 — 

Substitute these formulas into (fTTD : 

—27 sin 12p + 174 sin lOp — 420 sin 8 p + 378 sin 6 p + 153 sin 4p — 468 sin 2p 
° “ 2048p9 ’ 

—27 cos 12p + 84 cos lOp + 106 cos 8 p — 764 cos 6 p + 1099 cos 4p — 344 cos 2p — 154 


Bn = 


Cn = 


2048p8 

—27 sin 12p + 48 sin lOp + 140 sin 8 p — 336 sin 6 p — 71 sin 4p + 512 sin 2p 


1024p7 


Calculate the discriminant of equation ([ 12 ]): 


Do = Bl- AAqCo = (6561 cos 24p - 52488 cos 22p + 128628 cos 20p + 83592 cos 18p 
- 987134 cos 16p + 1543976 cos 14p + 702372 cos 12p - 4646312 cos lOp + 3755087 cos 8 p 
+ 3053616 cos 6 p - 4805144 cos4p - 4176688 cos2p + 5393934)/(8388608p^®). 


We used wxMaxima 12.04.0 for calculations. 

Obviously, Ag(A) 7 ^ 0, Dq{\) 7 ^ 0, so according to Lemma [2l the roots of equation ffTOj) in 
the general case have different asymptotics: 


M 2 '(A) 


sin p 


M 2 (A) 


1 + 6 cos^ p 
3 sin p cos p 


Since M 2 (A) = , only the root M 2 (A) is the required one. 

Finally, one can easily hnd M^{\) and solve classical Sturm-Liouville inverse problems by 
Weyl functions on the edges 62 and 65 . 

Now let us consider the case when the potential is known a priori on two edges. If they are 
Cl and 64 , then only two spectra Ag and A2 are sufficient to recover the potential on the whole 
graph. Indeed, one can solve IP(2), then apply Theorem [T] to the vertex V 3 , hnd gs and then 
similarly hnd q^. However, the knowledge of qi and ^2 do not allow us to recover the potential 
from two spectra by our method. If we have only Aq and A4, we can not recover q^. Similarly, if 
we know q^ initially, the knowledge of the potential on one of the boundary edges do not allow 
us to reduce the number of given spectra. Thus, if the potential is known on multiple edges, 
the number of required spectra depends on the location of these edges. 
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